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Abstract 

In the present paper, an empirical study of LIBOR (London Interbank Offered Rate) 
data is presented. In particular, a data set of interest rates from 1997 to 1999, for 
two different currencies and various maturities, is analyzed. It turns out that the 
random behavior of the daily increments for the interest rates series is non-Gaussian 
and follows a leptokurtic distribution. 
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1 Introduction 



In financial theory and practice, interest rates are a very important subject 
which can be approached from several different perspectives. 

The classical theoretical approach models the term structure of interest rates 
using stochastic processes and the no arbitrage argument. Various models have 
been proposed and can be found in Pagan (1997), Bouchaud (1997), Rebon- 
ato (1998), Wilmott (1999), Hull (2000), Mantegna (2000), Dacorogna (2001), 
Vasicek (1977), Cox (1985), Heath (1992), Hull (1993) [1-11]. Although they 
provide analytical formulas for the pricing of interest rate derivatives, the im- 
plied deformations of the term structure have a Brownian motion component 
and are often rejected by empirical data (see Chan (1992) [12]). The inadequa- 
cies of the Gaussian model for the description of financial time series has been 
reported since a long time ago by Mandelbrot (1963), sec also Mandelbrot 
(1997) [13,14], but thanks to the availability of large sets of financial data, 
the interest on this point has risen recently in Mantegna (2000), Dacorogna 
(2001), Sornette (2000), MuUer (1995), Mantegna (1994), Ghashghaie (1996), 
Mantegna (1997) [6,7,15-19]. In particular, the fat-tail property of the empir- 
ical distribution of price changes has been widely documented and is a cru- 
cial feature for monitoring the extreme risks (Embrechts (1997), Embrechts 
(1999), McNeil (1998), McNeil (1999) [20-23]). Many of the recent studies con- 
cern high frequency data of stock indexes or exchange rates (MuUer (1995), 
Ghashghaie (1996), Mantegna (1997), MuUer (1996) [16,18,19,24]). 

Also for interest rates, up to now, no universally accepted theory has been 
obtained for the description of experimental data. Dacorogna has investigated 
the intraday behavior of the interest rate markets (MuUer (1996), Piccinato 
(1997), Ballocchi (1997), Ballocchi (1999) [24-27]). Bouchaud et al. (Bouchaud 
(1997), Matacz (1999), Matacz (2000) [28-30]) have considered the interest 
rate curve as a vibrating string subject to random shocks along its profile. 
Along these lines, Sornette (1998) [31], has derived a general condition that 
the partial differential equations governing the motion of such string must 
obey in order to account for the condition of absence of arbitrage. Jean Nuyts 
(1999) [32] has developed a phenomenological approach based on Pade Ap- 
proximants. In the work of Jamshidian (1997) [33], a self-contained theory is 
presented for pricing and hedging LIBOR and swap derivatives by arbitrage. 
Appropriate payoff homogeneity and measurability conditions are identified 
which guarantee that a given payoff can be attained by a self-financing trad- 
ing strategy. An important recent development in the pricing of interest rate 
derivatives is the emergence of models that incorporate lognormal volatilities 
for forward LIBOR or forward swap rates while keeping interest rates stable. 
Zhao et al. (1999) [34] have introduced methods for discretizing these mod- 
els giving particular attention to precluding arbitrage among bonds and to 
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keeping interest rates positive even after discretization. 

In this framework, we have empirically studied the probability density distri- 
bution of LIBOR, in order to characterize the stochastic behavior of the daily 
fluctuations. 

In Section 2, we present the data set. The data are then analyzed and the 
results are given in Section 3. Section 4 contains a discussion and some final 
remarks. 



2 The data set 

LIBOR stands for the London Interbank Offered Rate and is the rate of in- 
terest at which banks are willing to offer deposits to other prime banks, in 
marketable size, in the London interbank market. 

BBA (British Bankers' Association) LIBOR [35] is the most widely used 
benchmark or reference rate. It is defined as the interest rate at which de- 
posits are perceived to be generally available in the London interbank at 11 AM 
(Greenwich time). 

Moreover, it is used as the basis for settlement of interest rate contracts on 
many of the world's major future and option exchanges as well as most Over 
the Counter (OTC) and lending transactions. 

The BBA, advised by senior market practitioners, maintains a reference panel 
of 16 banks. The purpose is to produce a reference panel of banks reflecting 
the balance of the market by country and by type of institution. Individual 
banks are selected within this principle on the basis of reputation, scale of 
market activity and perceived expertise in the currency concerned. The BBA 
surveys the panel's market activity and publishes their market quotes on- 
screen. The highest four and lowest four market quotes are disregarded and 
the remaining eight quotes are averaged: the resulting spot fixing is the BBA 
LIBOR rate. The quotes from all 16 panel banks are pubhshed on-screen 
to ensure transparency. BBA LIBOR is compiled each working day by an 
electronic vendor (currently Dow Jones Teleratc) and broadcast through ten 
international distribution networks including Reuters, Dow- Jones Telerate, 
Knight- Ridder, and Bloomberg. 

BBA LIBOR fixings are provided in seven international currencies: Pound 
Sterling, US Dollar, Japanese Yen, Swiss Franc, Canadian Dollar, Australian 
Dollar, EURO. LIBOR rates are fixed for each currency at monthly maturities 
from one month to 12 months. Rates shall be contributed in decimal to at least 
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two decimal places but no more than five. 



In the following, we have analyzed a data set of LIBOR interest rates r{T, t), 
where T is the maturity date and t the current date, for EURO and Pound 
Sterling. These data are shown in Fig.l and Fig. 2 where t goes from January 
2, 1997 to September 17, 1999, and T assumes the following values: 1, 3, 6, 9 
and 12 months for the Pound Sterling and 1, 3, 6, 12 months for the EURO. 

In Figs. 3, 4 the 1-month LIBOR is compared to the interest rates fixed by 
Central Banks at that time, namely the REPO (repurchase agreement) data. 
BBA LIBOR follows the trend determined by the decisions of Central Banks. 
In order to roughly eliminate these trends, in Figs. 5 and 6 the interest rates 
differences Ar{T,t)=r{T,t + At) - r{T,t), with At being 1 day and r=l 
month, are plotted as a function of the current date, for the EURO and the 
PoTind Sterling, respectively. A similar behavior is also found for the other 
maturities. Some large oscillations of Ar are induced by Central Banks. In 
any case, Ar heavily fluctuates around zero. 



3 Results 

In this section, we focus the attention on the probability distribution behavior 
of the interest rates increments Ar{T,t) defined in the previous section. To 
this purpose, we estimate ^(Ar), the complementary cumulative distribution 
function of the daily interest rates increments, defined as: 



where p is the probabihty of Ar(T, t). 

Because LIBOR data are supplied with only few decimal digits, it is interest- 
ing to examine the effects of different data cut-offs in the behavior of ^'(Ar). 
In Fig. 7, we plot the complementary cumulative distribution function for a 
simulated Gaussian stochastic process using data characterized by three dif- 
ferent decimal digit precisions. It turns out that the numerical rounding does 
not influence the results. 

Figs. 8 and 9 show the tail distribution behaviors in the case of EURO and 
Sterling Pound, respectively. In particular, in Fig.8, we report the empirical re- 
sults obtained estimating the probability density function of both positive and 
negative LIBOR increments with At=l day and T=l month. These empirical 
curves are slightly asymmetric and the negative variations are more probable 
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than the positive one. In the same figure, these two curves arc compared with 
the equivalent (i.e. with the same average and standard deviation) Gaussian 
complementary cumulative distribution. The non-Gaussian behavior is also 
evident from Fig. 9. In both Figs. 8 and 9, the empirical LIBOR data exhibit 
a fat tail or leptokurtic character, which is present for the other maturities 
as well. These observations indicate that the random behavior of Ar(T, t) is 
non-Gaussian and that using a Gaussian probability density function leads to 
underestimating the probability of large fluctuations. 

For a better understanding of what kind of stochastic process we are dealing 
with, we present some results on the power spectral density behavior (Kay 
(1981) [36]). 

The power spectra, S{f), for both r(l,f) and Ar(l,t) are reported in Figs. 10 
and 11 for the EURO and in Figs. 12 and 13 for the Sterling Pound. For r the 
spectral density shows a power law behavior. A linear fit gives a slope value 
a = -1.80 ± 0.02 and a = -1.79 ± 0.01 for EURO and Sterling Pound, re- 
spectively. A similar result holds for the other maturities. Therefore, we argue 
that the power spectrum analysis for r(T, t) indicates a stochastic process with 
spectral components decreasing as S{f) ~ (Feller (1971) [37]). As for the 
increments, in Figs. 11 and 13, the power spectrum is fiat, typical of a white 
noise process without seasonality peaks. These results are also corroborated 
by a similar analysis performed on Eurodollars interest rates for a longer time 
period (Di Matteo (2001) [38]). 



4 Discussion 

In the previous section, we have shown that the daily increment of the interest 
rate series is non-Gaussian and follows a leptokurtic distribution. The problem 
arises of how derivatives written on interest rates can be evaluated, given that 
the usual Gaussian white-noise assumption of many models is not satisfied. 
Indeed, a first partial answer, is that the central limit theorem ensures that, 
after a sufficiently long time, the increment distribution will tend to a Gaus- 
sian distribution. However, if the time horizon of derivative evaluation is not 
appropriate, deviations from the Gaussian behavior may lead to a dramatic 
underestimate of large increments with a consequent improper risk coverage. 

Although well studied in mainstream finance, this problem has received much 
attention in recent times, within the community of physicists working on fi- 
nancial problems (Bouchaud (1997), Bouchaud (1994) [2,39]). In particular, 
Bouchaud and Sorncttc (Bouchaud (1994) [39]) suggested the direct use of the 
historical probability measure, rather than the equivalent martingale measure 
for evaluating options. In this way, one gets an option price depending on 
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the expected rate of returns, a conscqTicnce which is not fully desirable due 
to the STibjcctive character of that rate. Assessing trends is a difficult task, 
as they depend on decision taken by Central Banks (as shown in Figs. 3 and 
4) and are based on macroeconomic effects. Thus, martingale methods could 
prove more reliable. As early as 1977, some years after the seminal paper 
of Black and Scholes, Parkinson generalized their approach to option pricing 
and explicitly took into account leptokurtic distributions (Parkinson (1977) 
[40]). More recently, Boyarchenko and Levendorskii have studied the problem 
of option pricing in the presence of a specific distribution which seems to fit 
well the empirical data in many instances: the truncated Levy distribution 
(Boyarchenko (2000), Koponen (1995) [41.42]). Along these lines, we believe, 
it will be possible to develop a consistent option pricing model taking into 
account the leptokurtic character of the empirical short to mid-term interest 
rate distributions. 
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Fig. 1. LIBOR interest rates r{T,t) as a function of the current date t, for the 
maturities T=l, 3, 6 and 12 months and EURO currency. 
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Fig. 2. LIBOR interest rates r{T,t) as a function of the current date t, for the 
maturities T=l, 3, 6, 9 and 12 months and Pound Sterhng currency. 
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Fig. 5. 1-month LIBOR increments as a function of the current date for EURO. 
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Fig. 9. The complementary cumulative distribution function of LIBOR increments 
for Sterling Pound. T=l month and At=l day. 
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Fig. 11. Power spectrum of Ar(T, for the EURO with T = 1 month and At=l 
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Fig. 12. Power spectrum of r(T, for the Sterhng Pound with T = 1 month. 
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